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ABSTRACT 

The tree order power spectra of primordial inflation depend upon the norm- 
squared of mode functions which oscillate for early times and then freeze 
in to constant values. We derive simple differential equations for the power 
spectra, that avoid the need to numerically simulate the physically irrelevant 
phases of the mode functions. We also derive asymptotic expansions which 
should be valid until a few e-foldings before first horizon crossing, thereby 
avoiding the need to evolve mode functions from the ultraviolet over long 
periods of inflation. 
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1 Introduction 



Cosmological perturbations from primordial inflation [TJ [2j [3] have a crucial 
importance for fundamental theory because they are first quantum gravita- 
tional effect ever detected [U [5], because they provide strong evidence for 
primordial inflation [HI [TJ [8] , and because they set the initial conditions for 
structure formation in cosmology [91 [101 [11] . Two sorts of perturbations are 
produced during single-scalar inflation: a scalar perturbation characterized 
by the field ((t, x), and a tensor perturbation characterized by the transverse- 
traceless field hij(t,x). The scalar signal has been imaged in the anisotropies 
of the cosmic microwave background [121 [T3], and by measuring the matter 
power spectrum with large scale structure surveys [2]. The tensor signal has 
not been imaged so far [121 [13] but strenuous efforts are underway to detect it 
through the polarization of the cosmic microwave background [151 US El HE] • 
The two perturbation fields are defined on the homogeneous, isotropic 
and spatially flat geometry characterized by scale factor a(t), with Hubble 
parameter H(t) and slow roll parameters e(t) and r)(t), 



ds 2 



(2 H 

-dt 2 + a 2 (t)dx-dx , H(t) = - , e(t) = , r)(t) = e - 
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Only the quadratic parts of their Lagrangians are relevant to current mea- 
surements (see [H] for an explanation of the full formalism), 
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The reported results [121 [131 E! f° r the scalar and tensor power spectra are 
consistent with evaluating the following 2-point correlators long after the first 
horizon crossing time tk such that k = H(tk) o(tjfc), 
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At tree order the correlators (SHSJ) can be expressed in terms of the scalar 
mode function v(t, k) and its tensor cousin u(t, k), 



A^(t,A0 = f-x4nGx\v(t,k)\ 2 + 0(G 2 ), (6) 
A 2 h (t,k) = ^x 327rGx2x\u(t,k)\ 2 + 0(G 2 ) . (7) 
The relevant mode equations and normalization conditions are, 

i 2 

v+ (3H + -)v + —v = , W v = vv* — vv* = — ^-r , (8) 

A; 2 i 
u + + — u = , = uu* — iiu* = — . (9) 

a 2 a 6 

The initial conditions for the mode functions (which correspond to the choice 
of Bunch-Davies vacuum) derive from the WKB solutions in the far ultravi- 
olet (that is, for k > H(t)a(t)), 

expl — ik f* dt' I 'a(t')] \ ^/Ha\ i 
v{t,k) — > , J / = J-^lJl + 0( — U, 10 

V ^ v/2A;a 2 (t)e(t) \ \ k J 

u(t,k) expHy^tVa^f /ffax. 

^/2kaHt) \ \ k ) ( 

As is evident from these asymptotic forms, both mode functions oscillate 
and fall off in the far ultraviolet. They become approximately constant near 
the time of first horizon crossing (One can infer the existence of constant 
solutions quite generally when the ^f^y terms in equations (JEHU]) become irrele- 
vant.) It is these constant amplitudes which determine the crucial theoretical 
predictions for the power spectra through relations (EE)- 

One of the frustrating things about this formalism is the need to employ 
approximation techniques, even for evaluating tree order formulae such as 
(JBHZj)) because equations (JHJ|9]) for the mode functions cannot be solved ana- 
lytically for general scale factor a(t) [20||2~T]. Examples of such approximation 
techniques include 



Assuming the slow roll parameter e(t) is constant; 

Matching the leading ultraviolet and infrared forms at t — and 
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• Employing the full WKB solutions. 

When greater accuracy is needed, one must resort to numerical evolution 
of ([HUH]) from the known initial conditions (llOHlip until well past the time 
of first horizon crossing [23]. Excellent numerical solution techniques exist 
[25j |26| 127] , but a large fraction of their power is wasted on reproducing the 
oscillations of the mode functions, which contribute nothing to the power 
spectra. This can be especially time-consuming for models in which there 
is an extended phase of inflation, or when scanning properties of classes of 
models. 

A more effective technique was developed recently for the tensor power 
spectrum [28] in order to work out the gravitational wave signal from a novel 
model of inflation in which the Hubble parameter oscillates for a brief time 
[29l [30] . The key to the new technique is to convert relations (jHJ) into an 
equation for the norm-squared of the mode function, M(t,k) = \u(t, k)\ 2 . 
This is the quantity that actually enters the tensor power spectrum ([7J), 
and it has a much more sedate evolution than the mode function. Evolving 
M(t, k) avoids the wasted effort of numerically simulating the oscillations of 
u(t, k). Further, an excellent asymptotic expansion can be derived for M(t, k) 
which must be valid until only a few e-foldings before first horizon crossing, 
no matter how long inflation persists [28J. It is therefore only necessary to 
numerically evolve M(t, k) for a few e-foldings, from just before tf. until it 
becomes constant. 

The purpose of this paper is to extend this more effective technique to the 
scalar mode functions using a trick [21] for converting solutions to u(t, k) into 
solutions for v(t, k). To further simplify the formalism, we express our results 
in terms of the tree order power spectra (JMZJ) directly, without ever mention- 
ing the mode functions. Section [2] derives the key equations for A^(t, k) and 
A^(t, k), and section [3] gives their asymptotic expansions. Our discussion 
comprises section |U 

2 Differential Equations for A| and A 2 h 

We begin by reviewing how one derives an equation for M(t, k) = \u(t, k)\ 2 . 
The first step is to write out the first and second time derivatives, 



M(t,k) 
M(t,k) 



uu* + uu* , 

uu* + 2uu* + uu* . 
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One next eliminates u and u* in f)13p by using and then recognizing 
factors of M and M, 

2k 2 

M(t, k) = -3HM —M + 2uu . (14) 

a 2 

The final factor involving uii* can be expressed in terms of M by subtracting 
the square of the Wronskian W u in expression from the square of M in 
expression (1T2]) . 

1 M 2 1 

- ^ = 2M X 2 * = M2 + ? = 2M + 2M^ ■ < 15 > 

Substituting (fl5|) in relation (111)1 gives the desired equation for M(t, fe), 

9P /if 2 1 

M + 3 ™ + ^M = |L + ^. (16) 

Deriving an equation for N(t, k) = |v(t, A;)| 2 entails the same first step, 

N{t,k) = vv* + vv*, (17) 
N(t,k) = vv* + 2vv* + vv* . (18) 

We next employ expression (|Sj) to eliminate i; and v* in favor of N and TV, 

/ f\ 2k 2 
N = -[3H + -)N -N + 2vv*. (19) 



r ' a 2 



Subtracting the squares of the N from ( TTfl) and the Wronskian W„ in equation 



1 N 2 1 

W-W? = 2^x2™' = /^ + _ 2«* = • (20) 

And the final equation for N(t, k) comes from substituting (120)) in (fT9")) . 

J v + ( 3H + f)* + ^ = ^ + _^_. , 21) 

Because equations (fl6|) and (12T|) are nonlinear, they incorporate the 
Wronskian normalization as well as the evolution equations for the mode 
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functions. Converting them to equations for A^(t, k) and A\(t,k) requires 
only changing the coefficients of the final terms, 




Al + 3HAI + %A\ 



(A 2 n ) 2 2G 2 k & 1 
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(22) 



(23) 



A 2 7r 2 a 6 A 2 ' 



3 Asymptotic Expansions for A^ and A 2 h 

Equations (T^2j) - T23|) require initial values for the power spectra and their first 
time derivatives. Of course these derive from the ultraviolet limits fllOHlip . 
which imply the zeroth order terms for an expansion in powers of (Ha/k) 2 , 



We shall develop the expansion for A\{t, k) to second order, and then employ 
a trick for converting tensor mode functions to scalar ones [21] to obtain the 
corresponding expansion of A\{t,k). These expansions are so accurate that 
there seems little point to numerically evolving for anything except the last 
few e-foldings before first horizon crossing. 

From the asymptotic limit (1251) one can see that the "big" terms in the 
evolution equation (|2"3"1) are those with no derivatives. It is best to segregate 
these terms, 




(24) 
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The desired expansion is now easy to read off [28J, 




Note that this expansion is in powers of (-^f) 2 , rather than the expansion 
in powers of 4^ one gets for the mode function. The better convergence for 
small jr is one more indication of the superiority of our method. 

The trick for converting A^(t,k) into A^(t,k) involves simultaneously 
changing the scale factor and the meaning of time [21], 



a(t) 
d_ 
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y/e(f) x a(t) 
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One can see that these two replacements convert the u(t, k) mode relations 
([nD into the relations for v(t,k). Replacements (I28H29P also convert the 
A 2 (t, k) equation (1231) into expression (I2"2"j) for A^(t, k). We could have used 
them to derive that equation had the straightforward derivation not been so 
trivial. 

To convert the asymptotic expansion ( |27l) into an expansion for A^(t, k) 
we first note that the replacements ( 12811291) imply the following geometrical 
replacements, 
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After many tedious manipulations the result is, 
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Unless e(t) or r]{t) is changing rapidly on the Hubble scale, the fractional 
error in AM]) should be of order (^) 6 - Had we used this formula to estimate 
A^(t, k) for t as close as one e-folding before first horizon crossing, the frac- 
tional error would only be about 2.5 x 10~ 3 ; for two e-foldings that would fall 
to about 6.1 x 10 -6 ; and the fractional error would only be about 1.5 x 10 -8 
at three e-foldings before horizon crossing. This is comparable to what loop 
corrections might give [3 1 J . 



4 Discussion 



We have derived nonlinear, second order differential equations fl22|) and (1231) 
for the correlators A^(t, k) and A^(t, k) whose late time limits give the scalar 
and tensor power spectra at tree order. Numerically evolving these correla- 
tors is more economical than evolving the mode functions (JEHU]) because one 
no longer has to keep track of the phases which drop out of the power spectra. 
This economy shows up in the asymptotic expansions and (1271) pertinent 
to the ultraviolet regime of k ^> H(t) a(t). Our expansions for the correlators 
are in powers of the small parameter (4^) 2 , whereas the analogous expan- 
sions for the mode functions are only in powers of In fact our asymptotic 
expansions (I34H27P should be so accurate that they could be used to provide 
the initial value data for A^(t,k) and A\(t,k) only a few e-foldings before 
the time tk of first horizon crossing. 

We have employed "Hubble parametrization" , in which the scale factor 
a(t) is assumed known and results are expressed in terms of background 
geometrical quantities H(t), e(i) and rj(t) defined in expression (pQ). This is 
especially useful for the tensor power spectrum because gravitational waves 
depend only on the background geometry, no matter what caused it. Many 
people prefer "potential parametrization", in which the scalar potential V(<f>) 
is assumed known (with the field starting at rest from some fiducial value) and 
results are expressed as derivatives of the potential. Approximate conversion 
formulae are, 



V'\ 2 1 



H 2 ~ -irGV , e ~ — — ( — J ,77 



3 ' l6irG\VJ ' ' 167rG 



v \v 



■ (35) 



Finally, we note that the same technology can be applied to derive an 
equation for the power spectrum of any field whose plane wave mode func- 
tions obey the Mukhanov equation J2JJ. Of course the associated asymptotic 
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expansion is also easy to develop. Both the equation and its solution in the 
regime of k 3> H(t) a(t) can be obtained from the analogous tensor rela- 
tions, f )23|) and ff27j) . by changing the scale factor and the time as we did in 
equations ([281129]). 
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